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Resummations of soft gluon emissions play an important role in many applications of QCD, among them jet 
observables and small x saturation effects. Banfi, Marchesini, and Smye have derived an evolution equation for 
non-global jet observables that exhibits a remarkable analogy with the BK equation used in the small x context. 
Here, this analogy is used to generalize the former beyond the leading N c approximation. The result shows 
striking analogy with the JIMWLK equation describing the small x evolution of the color glass condensate. A 
Langevin description allows numerical implementation and provides clues for the formulation of closed forms for 
amplitudes at finite N c . The proof of the new equation is based on these amplitudes with ordered soft emission. 
It is fully independent of the derivation of the JIMWLK equation and thus sheds new light also on this topic. 



1 Introduction 

The starting point for the issues discussed and developed in this paper is a beautiful exposition by Banfi, 
Marchesini, and Smye (BMS) pQ of an evolution equation describing interjet observables in hard QCD pro- 
cesses in which one asks for the energy flow E out into a region away from all hard jets. With E out much 
smaller than the hard scale of the jet, QCD radiation into this region is small and these observables offer 
an ideal opportunity to study nonperturbative effects. The central result of is an evolution equation that 
resums, for large Q/E out and in the large N c limit, all leading terms arising from large angle soft emission. 
In the following, I will refer to this equation as the BMS equation. Subsequent analysis of the consequences 
implied by this equation presented in the same publication confirm features found using numerical methods 
by Dasgupta and Salam The equation itself is nonlinear in the part describing contributions arising 

inside the jet regions and linear outside. The contributions inside, quite surprisingly, have a striking anal- 
ogy with the Balitsky-Kovchegov (BK) equation which resums small x corrections in the nonlinear 
domain. There it is used to describe saturation phenomena that occur in DIS at small x, possibly off large 
hadronic targets. For small angle emission in a certain frame, the analogy becomes even more striking: the 
contributions of BMS inside the jet regions become identical to the BK equation. That such an analogy, 
if it exists at all, is restricted to the in-region, is not surprising: in the derivation of the BK equation one 
assumes, in fact requires, the presence of non-linear effects everywhere and any geometric exclusiveness of the 
type imposed by jets was simply never considered. What appeared to be more surprising is the relation of 
jet physics to BFKL dynamics [3IHj, which dominates the BK equation in the small density limit. This link 
has been established in [§] and indeed comes into play in the limit of small angle emission in a judiciously 
chosen frame. 

I will be concerned with a different question: the generalization of the BMS equation to finite -ZV C based on 
the fact that such is already available on the BK side of the above analogy. 

The BK equation is the limit of a more general equation, the JIMWLK (Jalilian-Marian Iancu McLerran 
Weigert Leonidov Kovner) equation |l()lllllfT21ll^lll4llTnill6lllTlllSj 1 . a functional equation equivalent to an 
infinite hierarchy of coupled evolution equations known as the Balitsky hierarchy @H3j ■ 

BK emerges from JIMWLK by means of a factorization assumption strongly linked with the large N c limit. 
As both BK and BMS equations use eikonalized soft emission techniques, both in a strongly ordered setting, 
it is only natural that there will be a host of technical similarities in both situations that may serve as a 

lr The order of names was chosen by Al Mueller to give rise to the acronym JIMWLK, pronounced "gym walk". 



vehicle to derive a consistent generalization of the BMS equation that goes beyond the large N c limit. This 
derivation will be the goal of this paper. 

In Sec. |21 1 will shortly present both the BMS and BK equations and highlight their structural analogy. This 
will provide me with a tentative identification of the transition probabilities used in the BMS equation with 
averages of "dipole type" operators akin to those appearing in the BK equation, albeit in a different kinematic 
domain and geometrical setting. 

Using this identification in Sec. |3| I will suggest a nonlinear evolution equation that generalizes the BMS 
equation in the same way as JIMWLK generalizes the BK equation. From a single equation for a single 
function one generalizes to a functional equation, equivalent to an infinite hierarchy of equations for an 
infinite tower of correlators. 

I will identify real-emission and virtual correction parts of the resulting functional Fokker-Planck equation 
and discuss the issues of real- virtual cancellation and highlight the fact that this is in direct correspondence 
with infrared finiteness. 

Sec.0]will provide the translation of the resulting functional Fokker-Planck equation to an equivalent Langevin 
description. This will allow to identify explicitly both real and virtual contributions to the amplitudes in a 
closed form. At the same time one gains a tool to study this evolution equation at finite N c numerically, 
in analogy to what has been done already in the JIMWLK case (T^|. It will also provide insights into the 
structures that will be needed to derive the underlying Fokker-Planck equation. 

The remainder of the paper is devoted to a derivation of this equation from first principles. Due to the 
functional nature of the equation, the most economical way to proceed is to consider generating functionals for 
the objects of interest. This will allow me to write closed functional expressions for all the relevant amplitudes 
and transition probabilities at finite N c in terms of a gluon emission operator whose form is motivated by 
the comparison with JIMWLK - the amplitudes found are of course the relevant QCD amplitudes of |27fl2~T] . 

In Sec. [S] I will define generating functionals for amplitudes that describe the transition of any bare jet 
configuration with an arbitrary but fixed number of hard particles that emit any number of soft gluons 
<7i . . . g n in the strongly ordered domain u>i 3> . . . 3> w n , all softer than the original set of hard particles. At 
this stage I will focus on real emission. My initial example will be based on a bare qq configuration. 

Sec. is used to formulate generating functionals for transition probabilities corresponding to the above 
amplitudes that will then serve as the objects for which the evolution equation is derived. 

This is done in Sec. by considering derivatives of the expressions for real emission with respect to the phase 
space boundary in complete analogy to the method of BMS. Virtual corrections are then included in the 
usual way with a coefficient fixed to ensure real-virtual cancellation in the infrared. The result is a proof of 
the equation suggested in Sec. Inland a verification of the fact that it indeed resums all finite N c corrections 
in the kinematic domain under consideration. 

Conclusions and a short discussion on possible cross fertilizations between the two fields will conclude this 
paper in Sec. |H| 



To present the analogy, let me give a short account of both equations, starting with BMS. I will closely follow 
the exposition in pQ. 

The observable considered by BMS is given by 



Here E out is the sum of all energies of soft gluons emitted into the region outside some predefined jet region, 
C ou t and da n the n-hadron distribution and ot the total cross section. The dependence on jet geometry and 
total energy Q is understood. The process is dominated by iterative soft emission from cikonalized gluons that 



2 The analogy of BMS and BK equations 
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are not deflected from their original trajectory. The natural kinematic variables are thus the directions of the 
emitters (a,b,k,p,p in what follows) and the corresponding energies. Phase space integrals correspondingly 
will involve integrals over solid angles and energies. 

A Mcllin transformation represents the cross section in terms of transition probabilities G ab {E, E ont ) 

/dv e uEmt 
— tt~- G pP (E,v~ ) ~ G pP (E, E out ) , (2) 

were a saddle point approximation has been taken. 

G a b{E, E ou t), or Gab(E) where the explicit E on t dependence is not needed, satisfies the BMS equation 

E8 E G ab (E) = J ^ a s w ab (k) [u(k) G ak (E) ■ G kb (E) - G ab {E)] , (3) 

or equivalently 

Ed E G ab (E) = -Ed E R { °\E) ■ G ab (E) 

^a s w ab (k)u(k) [G ak {E)-G kb (E)-G ab (E)\ . 

47T 

a s := a s N c /TV, Rf^(E) is called the single log Sudakov radiator for bremsstrahlung emission 

,(0W _ f E duj f d 2 n k 



R ab (E) = J —J -^±a s w ab (k)[l-u(k)] = A-r ah , (5) 
and 

(k) = ( PaPb ^ = 1 - cos 6 ab 

Wab[ ' ( Pa k)(kp b ) (l-COSflofcKl-COSflw.) U 

is the soft emission kernel. Conventions here are adapted to lightlike momenta and the energies, LO p ,u) q have 
been factored out in (pq) compared to the 4- vector product: p.q = u> p Lu q (pq) = Lu p Lu q (l — cos6 pq ). 

The jet geometry is encoded in the definition of u(k), which emerges from a Mellin factorization of the energy 
constraint in Eq. QJ: 

e(k ut - J2 u ) = j dV ^w ' u ( fc ) = 9i "( fc ) + e "^ out(fc) , (7) 

leCout mV i 

with the O functions having support inside and outside the jet regions respectively. In fact, BMS show that 
to good accuracy one may take u(k) to restrict the phase space integrals to the outside region for the Sudakov 
radiator term and to the inside region for the remainder of Eq. . 

For more details on the ingredients as well as the physics of this equation see mElElEnilS] an d references 
therein. 

At this point I want to highlight the strategy employed by BMS in the derivation of Eq. (0 . This evolution 
equation was derived from the knowledge of the structure of the real emission part alone. The contribution 
of the virtual corrections was included after the fact using the requirement of real-virtual cancellation in the 
infrared. 

The real emission part of G needed for this argument can be written as 

dujj cPfL 



GpT l) (E,E out ) = i+Y, U « 

n=l*' i=l ^ 



uij Air 



u{h)Q(E-Lu l )}W n (pk 1 ...k n p) (8) 



where the phase space of soft gluons is cut by E and W n is the large N c factorized version of the transition 
probability of a color singlet into (gg)i iar dg" oft . 

W n ( P k 1 ...k n p) = { f p) . (9) 

{pki){kik 2 ) ■ ■ ■ (k n p) 
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The line of argument then starts by taking a logarithmic derivative E3e of Eq. JSJ which uniquely fixes 
the quadratic term in Eq. This term is affected by the jet-geometry as signalled by the factor u(k). 
The virtual corrections are then taken into account by subtraction of a linear term with a coefficient that 
ensures infrared finitencss through real virtual cancellation. Virtual corrections occur globally, both inside 
and outside the jet region, thus the factor u(k) is absent. This yields the second term in Eq. (J3J). 

That this indeed ensures infrared finiteness is best understood after splitting off the Sudakov term as in 
Eq. Noting that the integration in the Sudakov term is restricted to the out region one recognizes that 
there is no danger of encountering any ill effects from the poles in the kernel in this term. For the other term 
one needs to prove (see £Q) that G aa {E) = 1. Then the quadratic and the linear terms cancel where the 
kernels diverge, rendering the expression finite. 

The BK equation on the other hand involves small x kinematics in which in the dipolc approximation is 
taken. The latter corresponds to the large N c limit of the above. While again one deals with hard leading 
particles that emit softer ones in an cikonal manner, they all share a common longitudinal direction, selected 
by the geometry of the collision event. Instead of directions away from the hard creation vertex of the jet(s), 
one is left with (2d) coordinates in the plane transverse to the collision axis that characterize the location of 
the hard emitters. To set the coordinate system in the small x case, I follow the common choice to represent 
the hard degrees of freedom as given by eikonal lines along the x~ direction (the direction of the probe). 
The target, which is separated from the probe by a large rapidity interval r = ln(l/x) (to leading accuracy 
it moves along the x + direction) then supplies the field these eikonal lines interact with. This plays a role 
similar to the soft emissions in the BMS case. Any x + dependence is frozen out by time dilation. The field 
variables then are given as 



CO 

:= Pcxpj -ig J dz~A+ H (x+ = 0,x,z~)} (10) 



where A^~ o{t (x + = 0, x, z ) is the soft/ Weizacker Williams field of the target. The basic object entering the 
DIS cross section is the dipole cross section, the impact parameter integral of 

N T;xy = (N xy ) T N xy = tr(l - UlU y )/N c . (II) 

N T;xy is the dipole function, N xy the dipolc operator and the average (. . .) T is to be thought to arise from 
an average of the [/-configurations characterizing the scattering with a given target at a given r. The r 
dependence is target independent and IR safe and hence perturbatively calculable in form of an RG equation 
known as the BK equation. As usual, the nonperturbative information is relegated to its initial condition. 
To expose the similarity to the BMS equation it is useful to write the BK equation in terms of 

S T ; X y (S xy ) T S xy = tx(UlU y )/N c (12) 

instead of N T . xy . Then it takes on the form 

d T S T ; Xy 2^.2 J d, Z K^ X z y {Sr-.xz^r-.zy S T - Xy ') (1^) 



where 



is the BK-kcrncl. 



- (x - yf 

K "* v - (x-z)2(z-y)* ( ' 



Ignoring for the moment the Sudakov radiator term in Eq. fljl, the similarity of Eqns. Q and (|13fl is striking. 
A one to one relationship of structures emerges if one maps directions onto transverse coordinates, the kernels 
onto each other and tries to view G a b{E) as the average of some tr(U^Ub)/N c (where a and b represent the 
directions of the leading hard particles): 

G ab (E) (tt(UlU b )/N c ) E . (15) 
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Note that this is fully in line with the requirement that G aa (E) = 1 needed to ensure real virtual cancellation 
in the second term of Eq. Q. While this appears to be quite intuitive and fully in line with the physical 
interpretation it is not clear how precisely to perform that latter part of the translation. One would ask, for 
instance, how in detail to arrive at a consistent definition of these Us and the averaging process suggested by 
Eq. i|15|) in the light of the tree like branching process underlying the physics of the BMS equation. Additional 
questions would be how to reconcile their role in probabilities with that in amplitudes and many more in the 
same vein. 

Despite these open questions, the analogy is too strong to ignore. That the interpretation proposed in Eq. (|15(l 
should be possible suggests itself even more strongly, if one uses the observation made by BMS that in the 
small angle emission limit in a carefully chosen frame where the measure turns flat and the kernels agree 
completely: 

w ab {k) -► Km . (16) 

This is the region where Marchesini and Mueller [S] have established the link between jet and BFKL-dynamics. 

The benefit of demonstrating such a correspondence should be obvious to any reader familiar with the color 
glass condensate (CGC): If there is such an interpretation of the objects entering the BMS equation, then 
there is hope one might generalize it along similar lines as the BK equation, which can be viewed as the 
factorized limit of a more general functional equation, the JIMWLK equation which I will turn to next, in 
order to show how it reduces to BK under certain assumptions. This will be the analogy used to suggest a 
generalization of the BMS equation. 



3 From JIMWLK to BK: generalizing BMS by analogy 

To write the JIMWLK equation, consider the averaging procedure entering the definition of S T - xy to be 
performed using a weight functional for eikonal lines U: 

(...) T = / D[U]...Z T [U] (17) 



where D[U] denotes a functional Haar measure, in keeping with the group valued nature of the field variables 
U. 

The JIMWLK equation describes r dependence of the distribution functional Z T [U] . It was first written as 
a functional Fokkcr-Planck with a very compact evolution operator in |lfij and reads 

d T Z T [U] = -H FP Z T [U] . (18) 

The Fokker-Planck Hamiltonian is defined as 

1 

'2 



Hfp :=i iVJx&iV^ (19a) 



xt y -=-^ ]<PzK, xzy [{l-UlU z ){l-UtU y )} a0 (19b) 

and can be shown to be positive definite. In Eqns. (|19f) . both a summation and an integration convention is 
applied to repeated indices and coordinates. Where 



_ (a - z) ■ (z - y) 
xzv ~ (x-zf{z-yY 



is the JIMWLK kernel and is a functional generalization of the left invariant vector field on the group 
manifold [There is no coordinate derivative in the above] . Its precise definition will be given shortly. Let me 
first provide two structural comments: 

• Despite the apparent poles of the kernel at z = x, y, the Hamiltonian will lead to explicitly finite 
evolution equations because of the structure of the t/-dependcnce which vanishes at the potentially 
singular point. This is also the underlying reason for the cancellation of potential divergences in the 
BK limit. 
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• The equation is an equation for a distribution functional and thus a very compact way to represent an 
infinite number of (coupled) equations for individual correlators known as the Balitsky hierarchy 0]. 
To recover members of this hierarchy, say the evolution equation for a correlator 0[f7], one simply 
multiplies both sides of i|18|) with 0[U] and the averages with (|17|) . The result then is an evolution 
equation for (0[U]) T 

d T (0[U]) T = -(H FP 0[U]) T . (21) 

Because of the nonlinear nature of iJpp the r.h.s. will not be expressible by (0[U]) T alone but instead 
involves other correlators as well. The evolution equation of these new quantities will then also be 
needed, and the argument repeats itself, ultimctely leading to an infinite coupled hierarchy. 

Returning to the left invariant vector fields I note that for present purposes they may be taken to be defined 
as a variational derivative according to 

™% ■= -M/- ■ (22) 

su s is the ordinary functional (or variational) derivative w.r.t. the components of the U field: 

S 



U vM = SikSrfW (23) 



where Sxy '■= 5^ 2 \x — y) for compactness. Operationally this then leads to 

iVlUy := - U.rsW , iV a x Ul := PU^M . (24a) 

There is, of course a corresponding definition for the right invariant vector fields iVSJ,: 

iV%U v :=t a U x 8^ , iVJtft := -L^W . (24b) 

Their main properties are the commutation relations (which I display leaving the functional nature aside for 
a second) 

[iV a , i\7 b ] = if abc i\7 c [iV a , iV b ] = if abc iV c [iV a , iV b ] =0 . (25a) 

These then generate right and left translations respectively: 2 

e -w a (iv a )jj = Ue i" a t a e - iuja ^ a) U = e~ iu °' ta U . (25b) 

Functional forms of this of course involve an integral in the exponent: 

e -iL<(i^) Uy = U y e i < ta e"^" ( "->l/ 9 = r'<'°[/ s . (25c) 
V and V are interrelated by 

iV£ = -[UlViVi ; iVt = -[&]l b iVi (25d) 

and "representation conscious:" With the above definitions for the action on U and in the q and q 
representation, it automatically follows from representation theory that acting on U or U* in an arbitrary 
representation produces analogous formulae with the generators appearing on the r.h.s. in that representation. 

Where there are different U fields to distinguish, I will write instead of V£ for definiteness. 

These definitions and properties will become the core technical tool in the generalization of the BMS equation 
that is to follow below. For the time being, let me note that using both of these objects, £Tfp is elegantly 
written as 

Hfp = - ~ IC xzy [tV^iV; + iV a x rt a y + Uf(iV a x iV b y + *V^V b )] (26) 



Note that nevertheless, both correspond to what is called a left action of the group, as /(e I1J e lrl ,U) = f(e laJ ,/(e l71 ,U)) 
in both cases. 



G 



[integration convention for x,z,y]. While the factorized form of Eq. I|19|) is most useful in a derivation of 
a Langevin description of the evolution that allows a numerical implementation, this second form is more 
economical in the derivation of evolution equations for given correlators that follow as a consequence from 
Eq. JTHJl. 

Among these, the equation for the two point operator S xy is what I will turn to next in order to make contact 
with the BK equation. Multiplying both sides of Eq. I|18|) with S xy and taking the average over U according 
to Eq. (|1 T|> . one immediately arrives at 



~2tt 2 \ 

Us 

: 2tt 2 



tr(U m Ujh 
N c K 



(2fc X zy ^xzx fcyzy) ( ^Cj 



tr{U x Ul) 



N c 

ab tr{t a U x t b U y ) 



(27) 



In these expressions I have taken pains to label the contributions from the terms with and without an 
additional factor of U x by W and ' 2 ' for ease of reference. In the JIMWLK case there is no distinction 
between the kernels in the two cases and one finds that the linear combinations of JIMWLK kernels just 
assemble into the BK kernel: 



K, 



xzy 



(2K, (i) ~1C (i) -K (i) 

V ' xzy '^xzx yz i 



in both terms. To compare with the factorizing limit, one applies the Fierz identity 



[U z ] ab tr(t a U x t b uD = \(t?(U x Ut)tr(UzUl) ^tr^t/t) 



to rewrite Eq. QZjl as 



$r {,$xy)r 



a s N c 
2tt 2 



d Z K^xZy {S XZ S Z y S X y) T 



(28) 
(29) 
(30) 



This immediately reduces to Eq. (|30|) if one assumes factorization in the nonlinearity on the r.h.s. by replacing 



{Sxz $zy 



(31) 



Without this assumption, the evolution of the two point function contains a three point function, which in 
turn will couple to yet higher orders via an infinite hierarchy fully contained in Eq. (|18(l . The reason for this 
is the presence of U z in Eq. Q26[l. It is important to note that the factorization assumption (|31|) decouples 
and factorizes the whole hierarchy. Note that the real emission term originates solely from the contributions 
labelled (2) in Eq. I|27|l while the virtual corrections receive contributions from both. 

Since it would appear that the factorization assumption is closely linked with the large N c limit, it would 
appear that if an interpretation of G a b(E) in terms of averages over eikonal factors in analogy with S T . xy is 
possible, these objects should satisfy an evolution equation that is in close analogy to the JIMWLK equation. 
This is the premise on which I base the conjecture for the finite N c generalization of the BMS equation. To 
write it down I need 



1. to use the correspondences of variables and kernels noted in the above. In this step it proves useful, if 
I slightly redefine the functional aspect of the invariant vector fields to produce S functions adapted to 
the structure imposed by the solid angle measure. To this end I will write 

iV a p U p = t a U p 5(p - q) (32) 
where 5(p — q) contains Jacobian factors such that 

/ ^*(P "?)/(«) = Z(P)- (33) 
This will simplify expressions considerably. 
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2. More importantly, I need to take into account the inside-outside distinction of the jet geometry by 
allowing the counterparts of JC^ to be different for i = 1,2. Matching can then be done on the BMS 
level by considering the evolution equation for the two point function suggested in l|15|) . This is easily 
done and amounts to an algebraic exercise that fully determines both kernels from a rederivation of the 
BMS equation from its generalization via exactly the same steps displayed above for the JIMWLK/BK 
pair. Clearly both analogues of /C« will be proportional to w pq (k) with the proportionality factor 
f( l '(k) carrying the jet geometry via some uik) dependence. 3 

The result of this exercise is a Fokker-Planck Hamiltonian I will dub H ag (for "non-global"). I will present 
it in two forms, to parallel the two versions of the BMS equation, © and (0J. To this end I will introduce 
an additional function f^ 2 '(k) that will serve as the coefficient of the generalization of the Sudakov radiator 
term. One has the following forms of the Hamiltonian: 

H ng :=- g w uu (fc)[/ (1) (fc)(iV«iV$ + iV°iVS) +f (2 \k)[U k ] ab ( t V a ui \7 b v +iV a v iV b u )] (34a) 



in analogy with (J3J), and 

w m (k) [/« (k) (W» iV£ + iVJiVS 



H ng 



Us 

2tt 

+ f {2) (k){i^ a u iV a v + iV a JV a v + [U k ] ab (iV a u iV b v + iV a v iV b u ))] (34b) 

to parallel In both cases u, v and k are integrated over according to an "integration convention" that uses 
as its measure for any of the momenta. Retracing the steps leading from the definition of the JIMWLK 
equation to (|3L)fl by eye should make it obvious that indeed, the separation of terms in Eq. (|34b(l is such that 
the first line generates the Sudakov radiator and the second the nonlinear evolution inside the jet cones. 

A comparison of the result with BMS allows one to determine only the leading N c part of the /W. Already 
the 1/N C corrections are not controlled by the matching. To this accuracy the / are N c independent: 

f 2 \k) =u(k) (35a) 
/«(fc) =1 (35b) 
=(!-«(*)) . (35c) 

Actually, it would be unnatural for the Fokker-Planck Hamiltonian to contain any explicit N c dependence. If 
this were the case, it would require very special circumstances for it to incorporate all finite N c corrections to 
BMS in a consistent manner. For this reason I take Eqns. I|34|) . (|35t together with a Fokker-Planck equation 
of the form 

Ed E Z E [U} = -H ne Z E [U} (36) 

to define my conjecture for the finite N c generalization of the BMS equation. This result will, of course, be 
derived independently below. 

The most important feature of this equation is its infrared finitcness which directly relates to the real virtual 
cancellations being correctly encoded in Eq. I|34bjl . Considering the Sudakov term in H ng one finds that 
just as in the BMS limit, the phase space integration over k is restricted to the outside region where no 
hard particles are to be found. Therefore there is no divergence from the poles in w uv {k). The second term 
may be recast in analogy to Eq. (|19f) and the same finiteness argument used there clearly carries over to 
this situation. By looking at the BMS limit one sees that this same cancellation of terms corresponds to 
a real-virtual cancellation in the conventional sense. On the functional side, real and virtual contributions 
arise from different terms. Contributions with an additional factor Uk correspond to real emission, the rest 
to virtual corrections. 

The idea to interpret jet transition probabilities as correlators of path-ordered exponentials in singlet pro- 
jections properly takes into account gauge invariance, just the same way as the formulation of the JIMWLK 

3 Note that the calculations are somewhat simplified by the fact that w pp (k) = 0. The relation between kernels on the 
functional side to those on the BMS side as given by Eq. 1281 for JIMWLK/BK pair, becomes a matter of simple proportionality. 
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equation does. The equation suggested therefore would appear to consistently encode the physics it is in- 
tended to cover with a set of very strong constraints. From this perspective it would be rather surprising 
would it not be possible to arrive at this result from first principles. 



4 An equivalent Langevin description 

Translation of a Fokkcr-Planck equation to an equivalent Langevin description is a textbook topic, that for 
the present case is only complicated by the functional nature of the equation and the group valued nature 
of variables involved. That such a step is possible in the JIMWLK case has been first suggested in and 
worked out in detail in \22\ . 

Analytical results from the Fokker-Planck formulation arc generically hard to come by, although some infor- 
mation on the fixed point structure has been obtained in the case of the JIMWLK equation. This makes the 
alternative formulation even more useful, as it paves the way for a numerical implementation of the evolution 
equation. This has already been successfully carried out in the JIMWLK case in |19| . 

The Langevin formulation abandons the description in terms of weight functionals Z T in favor of one in terms 
of ensembles of fields which are governed by the corresponding Langevin equations. Explicitly, to calculate 
any observable O [U] of the fields U one writes 

(0[U]) T = Jd[U]0[U]Z t [U}<*± J2 °P\ ( 3? ) 

ueE[z T ] 

where, separately at each r, the sum is over an ensemble E[Z r ] of N configurations U. Its members can be 
thought of as created randomly according to the distribution Z T . Clearly, for N —> oo, the ensemble and Z T 
contain the same information. 

The dynamics are then encoded in a Langevin equation, a finite difference equation in evolution "time" r 
that contains a random force driven by a noise S. 

Prerequisites are (termwise) positive definitcness of the Fokkcr-Planck Hamiltonian in question. Only then 
can one introduce bounded noise integral - if needed separately for each positive term. If the Hamiltonian 
can be (termwise) separated into two conjugate factors in analogy to the structure of Eq. (|19|l . any correlation 
effects can be included in the Langevin equation. The noise can then be taken to be (termwise) uncorrelated. 
To establish the structure one may expect for the result, let me give the corresponding expressions for the 
JIMWLK equation. There, the Langevin equation schematically 4 reads 



d T [Ur^a = [U^n, [ J d 2 y E^ k [U T ]&* + <7« [U T ]\ (38) 



where 

is the "square root" of x, X% b y — f-xz^zy- Factorization is complete and no termwise split is necessary (or 
even possible without violating positivity). a is given by 

K\Ut\ := \%X% = i(~ Jd 2 zJC xz Mt a Ul, x U T , z )) . (40) 

The £ are independent Gaussian random variables with correlators determined according to 

(...> s = /£>[e](...)e-*«. (41) 



4 This is a continuous time version of the equation that strictly speaking is not unique. The path-integral derivation shown 
in 1221 and maybe more transparently in the appendix of 1191 makes it clear that we are to take a "retarded" prescription here 
in which the derivative on the l.h.s is taken as a finite difference and the fields on the r.h.s. are determined at the previous time 
step. 
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Fortunately all of the prerequisites listed above are also met for H ng . In order to achieve termwise factorization 
I first separate off the f( 2 ' terms and further treat the two contributions proportional to /W separately. All 
of these are positive due to the definition of u(k): 



u(k) = Q- m (k) + e-^Q out (k) 

ensures < u(k) < 1 and the same then for the coefficient functions in Eq. I|34b[) 
Correspondingly one would write the Hamiltonian on the JIMWLK level as 



(42) 



1 



where u, v are integrated over with ^Ll^Ll an( j 



ab 
Ami 



^-w uv (k) \fW(k)Q. + UtU v ) + f^(k)(l - UtU k )(l ~ Ulu v ) 

47T 7T L 



(43) 



(44) 



Factorization of x according to \ = is achieved by defining a three component structure representing 
the three separately positive terms alluded to before: 



pab;fi 



Correspondingly one has to introduce independent white noise for all of the components, 

s^ = {(e w )^,(€ (1 ' ) )^,« (2) )^}. 



(45) 



(46) 



The only thing left to cope with is the measure and the 5-functions in the noise correlator, which need to be 
such that 

(47) 



/ / f ^k f.nc;)i^c;/i\ / f pad;vr~d;i>\\ a 

\\ J l^^pk )\ J ~^ C ql "I )/-Xl» 



b 

m 



Now the measure is d co f 6d ^ and one needs the correlators to read 

= ^S(co S e p - co S e q )5(<b p - 0,)*V . (48) 
With these preparations and the identification r <-» In E one has a Langevin equation that reads 



k £;^[u T , p ] s*< 



47T 



(49) 



In comparison with the JIMWLK case, any sigma-terms vanish because w pq (k) satisfies w pp (k) = 0. 

Note again that any continuum notation is deceptive. The equation to consider is really a finite difference 
equation and upon iteration will allow for an interpretation of subsequent, ordered soft gluon emission. To 
this end it is instructive to rewrite the Langevin equation in terms of functional derivatives: 



9 T [U T - P ]ij — [U T .pit ]ij I £ pk ^[Ur-p] l3 t '.' 



Ir.fi 
k 



dfl„ dfl 



4ir 4ir 



is^[u r ., p ] st ; ; ^v^ ;g }[c/ r;; 



(50) 



The operator in the last version contains two terms in its third component, the component responsible for 
the in-region. Using a somewhat more compact notation for the t dependence they read: 



fdilg 


fdn k 


1 4tt j 


' "ItT 



dQ,„ dQ, k j a s p H 



4ir 4ir 



7r p.q 



a ; / / 



(51) 



Taking into account the fact that the Langevin description is defined only in a r-discretized sense, one might 
be tempted to conclude that the terms containing factors U q correspond to real gluon emission, while the 
others would generate virtual corrections, not only in this, but also the other components. 
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A formal solution to l|49|) is 



P1V 



P T exp 





fdflg 












1 4tt 


in 



c pq "fc <?J T , 



(52) 



which again is to be interpreted in a discrete sense. To write this expression I have again adapted the 
definition of the functional derivatives to include a (S-function in t: iV° p U T / q = —U T . p t a S(p — q)8(r — r'). 

This would have the interpretation of subsequent real emission with virtual corrections correctly taken into 
account. 

Instead of jumping to conclusions at this point, I will substantiate these ideas with a derivation from first 
principles that follows the strategy of BMS. For this one needs a formulation in terms of functionals capable of 
handling finite N c corrections. This I will turn to next. While doing so, a number of interesting structures, like 
a closed, finite N c version of Eq. (JHJ) will emerge alongside similar expressions for the underlying amplitudes. 



5 Amplitudes in the strongly ordered domain 

In order to find a suitable starting point to go beyond the large N c limit one has to go rather far back and 
start with a general discussion of soft gluon amplitudes along the lines given already in |20U21j . 

To help organize the argument, I will start with the definition of a generating functional for the (tree level) 
amplitudes of (real) soft gluon emission from a qq pair in the strongly ordered region 

^fc„ < ^fc„_! < • • • «w p < u q (53) 

written as 

n=0 J 

The amplitudes for n-gluon emission from a initial qq pair are denoted by -4Qp^'"fc" )j where the momenta 
q,p,k\ . . . k n and color indices i, j, a± . . . a n are explicitly listed, a corresponding set of Lorentz-indiccs /ii . . . fi n 
is suppressed. They are isolated by n,-fold (functional) differentiation w.r.t. H at S = 0. At this stage, 3 is 
just an external source, the relationship to the noise of the Langevin description will become clear later. 

The ^Qpfc 1 1 ' fc n ) ar e known to satisfy an iterative structure in which the color singlet — > qqg^ t ^amplitude 
follows by induction from the corresponding amplitude with n soft gluons. As a result, the 0-gluon term in 
Eq. (|54|) fully determines the whole functional, although the explicit construction of the n-th order expression 
with the tools available to date becomes more and more cumbersome with growing n. The rules of the game 
have been explicitly demonstrated in |21| and require no recourse to the 1/N C limit. 

I will now try to establish a method that will allow me to give a closed operator form of the full tower 
of amplitudes contained in Eq. I|54|l . While explicit generation of a term with arbitrary n will still be 
cumbersome, the fact that I have a closed mathematical expression instead of a prescription for an iterative 
procedure will allow me to use it in calculations and eventually in the derivation of the evolution equation. 

For this purpose, a slight generalization of the definition l|54f) is operationally somewhat more useful in that 
it allows me to specify an explicit functional form that casts the iteration step of |21j as a simple functional 
operation. This operation will involve nothing more complicated than the functional differentiation rules 
given in Eq. I|24() and a suitably defined n = term to start the process of generating the amplitudes. 

The starting point is the definition 

A« [U,E] :=E/^.- ^A^Z)^^ ■ ■ ■ UZ b ^t PtUUoh ■ (55) 

n=0 J 
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Here I have simply added a factor of U in the appropriate representation to each leg of the diagrams contained 
in the definition of the ordered amplitudes. Diagrammatically the first few contributions read: 






Uk 2 E k . 2 




+ . . 



(56) 



Obviously Eq. 1)54(1 is just the special case at U =1: A l J q [E] = A l J q [l, S]. The rationale behind this definition 
lies in the very nature of soft gluon emission underlying the construction of the amplitudes. Per definition, 
soft gluon emission does not change the direction of the parent (the emitter), in that role as a parent the 
entity (be it a quark or a gluon) is therefore written as a path ordered exponential along the direction of some 
momentum. Emission of a softer gluon, that in later steps will also serve as a parent will then necessarily 
add an adjoint factor Uk and an eikonal emission vertex Ji k t a (with J^ k := |^r), where the generator is to be 
taken in the representation of the emitting object. As all these objects are encoded as eikonal lines U, this 
will automatically occur if one writes it via invariant vector fields iV" just as in the rewrite of the Langevin 
equation at the end of Sec. 21 In functional form, one ends up with the single emission operator 



4-7T 47T 



(57) 



Note how closely this resembles the real emission part of the Langevin equation 1(49(1 . Eq. 1(55(1 should then be 
completely determined by the n = term in the sum which describes qq-p&ir without additional soft gluons: 
higher orders should just follow from repeated application of ((57(1 . This needs now to be checked against 
what is known about the soft gluon amplitudes. 

First one needs to give an initial condition for this iteration, the bare qq term. This reads 



[A<°>] 




:= M 2 (q,p)[UlU i 



(58) 



and corresponds to a zero order amplitude of the form A( l ^ p ) = M2(q,p)8 l: > 5 
The above prescription then claims that the qqg so ft term is given by 



[A«]« g [c/,H]={y. 



47T 47T 



9 J ik u k -fc * v ; 



(59) 



To demonstrate once how to use the machinery, I will go through the steps explicitly. First insert [A^J^f/, 5] 
from ||SSJ|. This yields 



[AW]« [U,E] = { / gJ? k U? ~ b k ^}M 2 { qiP )[UiU, 



q\n ■ 



(60) 



All that is left to do is to use the differentiation rules 1(24(1 (with the adaptation to phase space l(65|l taken 
into account) and differentiate the two factors U£ and U q . The variation eliminates the I integral and one is 

5 For M n I have adopted the notation of 121 11211 
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left with 



[A (1) ]^[C/,S] 





^M 2 (p, q )A; q {k)[U^U q \ l3 Ut^ 



where 

Apq{k) = J pk - J qk ■ 
Now one isolates the amplitude via a variation in 5 and reads off 

= M 2 (p, q )A^ pq (k)[n j =: M 3 (p, 




—k 



(61) 



(62) 



(63) 



with Ms(q,k,p) again as in the notation of [21)112 1| . Eq. H63|) in fact is identical to Eq. (40) of [2] and an 
additional gluon emission, i.e. iteration with Eq. I|57[) . leads to their Eq. (42) for the amplitude with 2 soft 
gluons: 



[A (2) ]^[[/,S] 












U kl 3 kl 




:i 3 kl + ~< 


\ ~ „ + ~^ 


U k2 E k2 








\u; 










u q 








Uk 2 "Bk 2 




Uk^kt + 


*\ - + ' 




Uk 2 ^k 2 




\\ Uk 2 ^k 2 


ul 


\ w 





(64) 



(The color algebra for the explicit comparison is given in Add. lAl case n = 1.) The pattern simply follows 
the rules of functional differentiation and a comparison of analytic expressions for the amplitudes as given 
by the iteration rules of |21j was performed explicitly using symbolic algebra tools up to n — 6, where the 
task starts to become time consuming. To compare beyond finite orders, a few additional tools are required. 

In the above I have implemented the ordering "by hand," simply by requesting the hierarchy (|53|l for the u ki . 
To write the all orders expression advertised above, I need one more ingredient: the addition of an energy 
integral to the emission operator (|57|> and an appropriate extension of the definition of iV p to iV° defined 
such that 

WZ p , p U q = t a U p ~d(p - q) (65) 

where 5(p — q) is adapted to the measure used in the momentum and energy integrations. The natural choice 
for the latter is to replace the solid angle integrations in i|57[) by full fledged phase space integrations: 



dflk 

47T 



d$ k := 



d^k 



6(k 2 )e{co k ) 



(2tt) 2 4tt 



0(Wfc) 



(66) 



S(p — q) is then defined as the appropriate S function on the forward light cone, i.e. with Jacobian factors 
such that 



dfaSip - k)f(k) = f( P ) . 



(67) 



(Note again the similarity of the roles played by r-ordcring in the Langcvin description and energy ordering 
in this case.) With these definitions it is straightforward to verify that a replacement of the solid angle 
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integrations in the above by phase space integrations leads to no change in the previous results, if one 
implements energy ordering for subsequent definitions with ordering 0-functions, just as one would for the 
contour parameter in case of path ordered exponential. To leading logarithmic accuracy this will not change 
the result for the evolution equation. 6 

By now it should appear to be the natural choice if I write for the generating functional (again suppressing 
the w-label for compactness) 

A* [tf, 5] =P Uk exp { J d^ k g J^tjf ^iVf}[A^ q [U, 3] (68) 

where P Wk implements the strong ordering of Eq. . I should actually use a notation that displays the de- 
pendence on the functional form supplied in the zero emission part and for instance write A^ q [U, a, M 2 (p, q)]. 

In analogy with path ordered exponentials, this implies that the exponential series is to be interpreted as 
P«i h exp { J d$ fc d$, g J? k Uf S'^iVf } 

-EMI f d **" d **» 9Jl\M: K ^ n ^t:}--{ [d* kl d$ h gJ&jjg* B^iV?}] (69) 

n=0 J J 

where the 1/nl are absent due to the explicit ordering of the subsequent soft emissions. To allow for more 
compact formulae below, I will introduce the notation 

U[H, U] := P Uk exp { j <f* fc d$, g J? k Uf S» ; "iVf } (70) 

for this functional differential operator. The explicit mention of the [/-field in the notation will be needed 
later when different types of cikonal fields will appear. U[S, U] will turn out to have a meaning independent 
of the bare (n = 0) term it is used to act on. It will be called the (real) shower operator as it will show up 
wherever one is forming that ordered soft gluon cloud around a hard seed. 

With these tools a direct comparison with the iterative scheme of j^J at all orders becomes possible. There 
it was shown (c.f. Eq. (45) of [SJ) that the n-th order amplitude is of the form 

KX'X)= E M n (q,k h ,...,k ln ,p)[t h &] tJ (71) 

n„ +2 (i) 

where n„ +2 (/) denotes the set of permutations of the {p, l\, . . . , l n , q}. Moreover, the amplitudes satisfy the 
recursion relation 

n 

A QZX a C + \) = J2 E M n (q,k ll ,...,k ln ,p)A ljh+1 (k n+1 )[^ ...t l H n+ H 1 ^ (72) 

j=i n„ +2 (i) 

(c.f. Eq. (48) of 123) ■ This determines all orders, once the first term is specified. 

Since agreement has already been established for the first few terms (n = 0,1,2 were compared explic- 
itly, terms including n = 6 have been checked), all that is left to show, is that the functional prescrip- 
tion given above leads to the recursion relation (|72|l . This is readily achieved: Assuming the n-th order 
term, [A^]* J g [J7, H] to have the form l)7ip. one simply uses a single soft gluon emission operator to find 
[AW]V[U,3\: 

d* kl ...d® kn A§x:li)K bl *l\ ■■■u k : K zt PlHU q ]h • (73) 



6 Also the results below for the structure of the multiple ordered soft emission can be obtained both ways - to have a 
formulation with all ingredients (including the ordering) written explicitly is merely a convenience from this perspective. What 
the formalism docs, is to allow to give a closed expression for the generating functional, Eq. 1681 . 
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It is straightforward to carry out the differentiation: The product rule leads to a sum over j similar to that 
in Eq. Tlien one uses [? a "+ 1 ] b "+ lQ " [t a "] k i = [t n+1 ,t a *"]ij to combine the J factors into the Ai jlj+1 (k n+1 ). 

After this rearrangement the sum only runs over j = 1, . . . , n instead of p, 1, . . . , n, q and matches Ij72(l . The 
algebraic details of this rearrangement are given in App. ^ This shows that one indeed creates the same 
amplitudes with the full color structure, i.e. without any recourse to the 1/N C approximation. 

While in the above, I have done nothing more but to provide a functional form of a (real) ordered soft 
emission amplitude, I will now show that there is structure contained in the above that is very similar to the 
information contained in the Langevin equation and will eventually lead to the RG equation stated earlier. 
The main observation in this regard is that, as a consequence of the strong ordering, any s Jh.in in which k is 
at (or near) the phase space boundary, there will be only a contribution from the first emitted gluon (n = 1). 
As a consequence, for such variations one gets 

U] = U[S, U]{ f J? k UfiVt} . (74) 

Note that 

Af qk [U, 3, M 3 (p, k, q)] := U[3, U] Uf[Ult a U q ] ij M 3 (p, k, q) (75) 

is the qqg counterpart to Eq. I|68JI . In short, Eq. (|74|l states, that under variation at the phase space boundary, 
one relates the generating functional for a given tower of amplitudes to that of a tower with an additional 
hard gluon. For the above example it relates the amplitudes with a hard qq to those with hard qqg content: 

-^A* [U, E, M 2 (p, q)] = tf b qk [U, 5, A£ q (k)M 2 (p, q)] , (76) 

or diagrammatically 




These two items, the ordered nature leading to (|74J) under "hard" variations, and the fact that in this situation 
the number of hard legs of the soft emission amplitudes is increased by one as exemplified in (|77(l , will be the 
core observations behind the evolution equation I am aiming at. 

6 Transition probabilities from S averages 

As a first step to understand how to translate the above into expressions for transition probabilities, let me 
take the situation of a (real) qqg" o{t soft emission amplitude in which one obtains the standard result 

<ilrZ)i A Cit:Z)V = E E M n ( q ,q h ,...,q ln ,p)M n (q,q h ,..., qin ,py 

n„ +2 (z) n„ +2 (Z') 

x tr(t a 'i t ai *t at 'n f l 'i) , (78) 

in which I have labelled the gluon momenta in the final state by ^i, . . . , q n . The same expression is obtained, 
if one considers the S averaged product of functionals 

([A<">]« [tf,3][A<">]« [U,E]Xs (79) 

if the average is defined with a Gaussian distribution for ^ in which only strongly ordered modes q\ , . . . , q n 
have support. I write 

(. . .)„ H := det(M„)5 / D[S] . . . e"^ d ^ d ^ s l M ^ 3 o (80) 
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where d(j>k denotes the phase space integral for momentum k and I have suppressed discrete indices. M 
defines the E correlator, which in the fully exclusive case considered above would read 

n 

= (3»'"H£") niB = S ab g^~S(p -q)J2 Kl ft) (81) 

»=i 

In this case, because of the strong ordering of momenta in the amplitudes, the E integral will strictly pair off 
the one available momentum fcj in the amplitudes falling into the same range as g,. 

As such, this is only true if I have ordering in the emission vertices, either done by hand in the 2-d version, or 
via Pu k in the 2+1-d formulation. Only then, unwanted cross terms are excluded. The U factors just cancel 
trivially. 

To make contact with the non-global observablcs of the BMS setting, one needs to depart from the fixed n 
situation and replace the VJ 5{q — q{), which selects a given set of final state momenta in the case where the 

i=l 

number of final state gluons is fixed, by a cutoff criterion that restricts real emission by geometry and energy. 
By writing 

Mpf b > v = (H^S^)e = S ab g^~6(p - q) d(E - u q )u(q) (82) 
one can indeed go beyond a final state with a fixed number of gluons and consider 

(A%[U,E]A% q [U,E]% , (83) 

which will provide the finite iV c -generalization of of BMS. 

In fact any degree of inclusivencss may be imposed by modifying the restrictions on the allowed modes 
appearing on the r.h.s. of Eq. Q82[l. 

At this point a few comments about the meaning of E as it appears in the above arc in order. Viewed from 
a diagrammatic point of view, the E stand for nothing else but the final state soft gluons. This explains why 
there is an additional factor of i together with the use of the free gluon phase space measure in the averaging, 
compared to structures encountered in Sec.0] There, E should be viewed simply as a perturbative gluon and 
the same ingredients as those encountered here would show up in the formulation of S'-matrix elements via 
the amputation of external legs. 

This realization should make it clear how to include virtual corrections right from the outset, by a fairly 
simple modification of the above. I will refrain from doing so here and instead keep a closer parallel with 
the strategy employed by Q][2n|[2T] and use the expressions for real emission to deduce the structure of the 
evolution equation. 



7 Evolution equations for soft semi- inclusive quantities 

For semi- inclusive quantities like the non-global observables of £Q, in which one sums over a given, limited 
phase space volume of soft gluons, it is natural to ask for the dependence on that phase space boundary. 
Unlike a direct calculation of the average in Eq. (|83fl this should lead to a tractable result that exhibits new 
structure. In fact, I shall demonstrate that one recovers the generalization to the BMS equation suggested 
above. 

To arrive at an evolution equation, I will simply take a derivative w.r.t. to the phase space boundary of the 
emitted gluons in the semi-inclusive probability l|83|l . This is in direct correspondence with the procedure 
used in BMS. Here it is essential that one includes the phase space boundary into the definition of M. 
Then the result is most transparently displayed using the following simple relationship, based on functional 
differentiation and Legendre transformation in the case of a Gaussian action ("free theory"): 

(W[E}) 3 =N I D[E] W[E] e-^r^ = W[-^j}e~^ ^ ^ = e~ ' ° »* ~° W[E ] ^ (84) 



j=o 



1=0 
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for any functional W[S]. Here notation has been condensed even further, with all momenta, integration 
signs and measures suppressed. 7 Legendre machinery is used for the last equality sign. The "classical field," 
So = iMJ : vanishes at J = 0. 

The canonical example for W[E] here of course is 

(W[E]) S - {A%[U,E]A% [U,E}% , (85) 

which, by construction, is the finite N c generalization of G^ al \ 

For this example one immediately realizes that the exponential of the second order differential operator on 
the right hand side will perform precisely the "sewing" implemented by the Gaussian weight in Eq. (|83|l as 
long as M is diagonal in energies as in Ij82(l . In this case, the energy ordering ensures that 

hh 1 ^^^^ = (4 A - [ ^ 3o] ) M (4 A - [f7 ' 5]t ) ■ (86) 

Let me now carry out the derivative with respect to the phase space boundary in analogy with the derivation 
of the BMS equation sketched in Sec. El With M = M(E) a function of the phase space boundary, one 
immediately finds a general expression for the logarithmic iS-dcrivativc of the above expectation value that 
reads 

Ed E (W[E])^E) = --e-^-s JL-Ed E M{E)^r W[E ] . (87) 

2 to^o to^o H =o 

The main point about this seemingly trivial exercise is its use in conjunction with what is already known 
about functional derivatives of the amplitudes of interest: Eqns. 174|) . 176|) . Note that 8eM{E) will force 
the variations in the factor taken down from the exponential to be at the phase space boundary. Choosing 
W[E) = A l J q [U,E]A l J q [U,E]\ the square of the generating functional of <?<7<?" f t amplitudes, this will lead 
to the appearance of the square of the generating functional of qqg <?™ oft on the right hand side as follows 
from (GU), E3. 

This marks an important difference compared to the BMS case at finite N c : the equation does not close. 
Instead higher and higher hard correlators will enter: Clearly, an evolution equation for the qqg <?™ oft ampli- 
tudes will then couple in turn to the corresponding object for qqg 2 g" oft amplitudes and so forth. One is faced 
with an infinite hierarchy. All quantities encountered in this hierarchy can be defined in complete analogy to 
the examples presented in H68JI and (|75|l . One simply uses the desired combination of [/-factors as the "bare" 
term in the tower of amplitudes one is interested in. 

This implies that firstly, it is not sufficient to consider only the evolution equation for qqg" o[t amplitudes 
alone. Instead one has to capture the complete infinite coupled hierarchy of equations in one. Secondly, it 
indicates that such a step is possible, since the objects of interest have already been identified and indeed 
can be written down in a general form. The evolution equations then follow from 18711 . 

I therefore proceed to define the (tree level) generating functional based on m hard particles, a general 
"antenna pattern," by 



A£;;%jU,E}Ai\:X[V,E? = U[E,U] (UV*)$ ® . . . ® (UV*)™ U[E,V? 



u=v 



The factor {UV^)p} ® . . . ® (UV^)^ in this expression generalizes the expression for the bare particles of 
the qq case (the A l J q [U, E)A l J q [U, S]t from above) from 2 to m hard legs, and the U[S, U] ... U[S, V}^ 
implement the showering. I have given only q and q factors as any gluon can be written in terms of these as 
(U~vrf k b = 2tiit a (UV^)lt\UV%). 

One can now study the evolution of the expectation value of a given object of the type (|88|l or directly study 
a generating functional for all such objects in one go. The latter is given by 

W[jt,j] :=(U[S,C7] eOtriW+oMt/vt)^ u[SjV] t^| ( 89) 



7 For compactness I will do so throughout this section. Where indices and momenta are shown it is with the understanding 
that there will be integration conventions for repeated momenta with d<E>fc as the measure. 
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with integrals in the source exponents understood. Appropriate variations with respect to j and j* at j = 
will then select a given tower of in — > hard<?™ oft probabilities entering the hierarchy. Since, besides the phase 
space constraints contained in M, it is the particle content of the bare terms that define the amplitudes, a 
notation that emphasizes this is needed. To this end I will write 



= :((^t)(t>...«([/yt ) (t) (90) 

where the UV^ average in the second line is to be taken with a yet unknown weight that reproduces the l.h.s.. 
The idea is to interpret the soft gluon average as an average over UV^ configurations. This is completely 
legal, if somewhat formal as long as the only definition of the weight of the average is through l|9()(l . 

Nevertheless, it is this reinterpretation that will allow to make contact with the Fokker-Planck formulation 
of Sec. |21 and that is all that is really needed. 

For m = 2 and an adapted choice for the color structure one encounters 

{ir{{UV%{UV%)/N c )^ =: G%*\E) (91) 

and thus provides the precise definition of what was still a bit vague in i|15|) • Eq. (|90|l is of course much more 
general than that: it defines the real emission part for a general antenna pattern, as these will all be needed 
for the evolution equations. Inclusion of virtual contributions will be discussed below. Note that the eikonal 
lines appearing in (|90J) have the interpretation of a product UV^ of contributions from both the amplitude 
and the conjugate amplitude. 

With JHlJl one then has 

W[i f ,j] ^e- ^ s ^^UfSo,;/] e i2trj ' t,7yt+l2tr(c/v ' +)t: '' U[E. ,V}^ = ( e i2trjWt-M2tr(Dvt)t^»^ _ 

(92) 

This demonstrates that ^N\j',j\ is in complete analogy with (the real emission part of) Z[J, jt] in the 
derivation of the JIMWLK equation |16j . 

To find the evolution equation for this (meta-) functional and with it the infinite hierarchy of evolution 
equations alluded to above, all that is left to do, is to put together Eqns. JH7J), l|86|) . and (|74|l to get the 
dependence on phase space boundaries for arbitrary [/^-correlators (bare m-jet probabilities). This yields 
the r.h.s. for the real emission contribution: 

.A. M{E) . ,.! 

e- — - — U[S ,£/] 

{g J&Ufi^yEdEMiE)}*^ [g J^ k ,V^\v a v , } e ^uv^ + ,2tr(uv^ 3 



w[jt, j] =(u[s, u] (c/vt)(t) . . . ® (uv^)M u[s, v]t 



u[s ,y]t 



(93) 

u=v 



To arrive at a meaningful answer one still has to supplement virtual corrections. They are correctly in- 
corporated as usual, by simply subtracting the corresponding term proportional to the original transition 
probability. This was done in the derivation of the BMS equation and it also applies to the full hierarchy. 
For concreteness, and without loss of generality, I will demonstrate how to achieve this for the case of m hard 
particles shown in Eq. I|9(J|) . 

To understand the constraints on the virtual corrections, let me first study the real emission contribution in 
some more detail. Omitting, for the moment, the operator that implements the soft gluon shower and the 
expectation value in the expression for the real emission part, 



■U[S ,C/]...U[S ,^] t 



u=v 
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one has to consider 

{9 J?M h Kl][Ed E M{E)] w 5 bb ' g^{g Jjf'^'iV^} (UV*)™ ® . .. ® (UV*)™ (94) 

and evaluate the variations. Inserting the BMS choice i|82[l for M, it is easy to collect all the factors and to 
make the k integration explicit (I keep the integration convention for I and I' for readability): 

m=J ^ESiE-w^^wuMiUV^in,^ (UV%»®...®(UV%» 

^^ Wll , {k)i uv\S^ b v , (c/yt)(t) 8 . . . {uv^I\^ =e . (95) 

In this expression u>i and uj[ will be made hard when the functional derivatives hit the U and V. The energy 
part of the corresponding phase space integrals will be eliminated at the same time. In summary one will 
end up with 3 "hard" momenta k, I, V . 

Carrying out the differentiations a sum over pairings will emerge, with different color structures depending 
if one hits q or q lines. Concentrating on the individual terms in this sum, it is obvious that one has to 
deal with 4 different combinations, corresponding to the pairings qq, qq, qq, and qq. This is a structure well 
known from the the derivation of both the Balitsky hierarchy |J and the JIMWLK equation in ^B] . There I 
had first to unearth the invariant vector fields from a diagrammatic calculation that has started off from an 
entirely different setup. Since this tool is already available, I will not go through the line of argument used 
there, but instead make use of the unifying properties of these operators right away without exploring these 
structures explicitly. 

This may be done by noting that Eq. (|94|l allows for a very elegant rewrite: 

m ~w pq (k)u(k){UVf)t b (^uv^\uvt )t + ™? uvt) W b {uvf)p ) (C/V/t)(t) ® . . . ® (UV%\1 . (96) 

Eq. (|96J) follows from a direct correspondence of ^ajv^) k w ^h i^u k acting on the amplitude and iV/^-ytv 
with «Vy- fc acting on the complex conjugate factor. In Eq. (|96|l it is possible to anticipate that the energies will 
all be made hard by functional differentiation on the hard factors (UV') Pi and thus to go back to the original 
definition of the variations without the energy 5-function. This allows to use an integration convention that 
employs the solid angle integrations only and a symmetric treatment of k,p,q in this equation. Reinstating 
the shower-operators and the expectation value, one arrives at the real emission part of the evolution equation 

(^p 9 (fc)^^(^t)f^V^ t)/ V^ vt)? +zV^ yt)9l V^ vt)p ) (UV*) p V®...®(UV^l) uvt (E) . (97) 

To discuss virtual corrections one needs to consider the physical process at hand. While the real emissions 
are confined to the inside regions and hence carry a factor u(k), virtual corrections appear everywhere and 
therefore have no such factor. Moreover the factor (UV')^ will be replaced by 5 ab and and the vector fields 
will have to act twice within either the amplitude or its complex conjugate, leading to i^ljjyi^ ^luv) anc ^ 
an analogous barred contribution instead of the mixed ones in Eq. I|97|l . For symmetry reasons one therefore 
expects the virtual corrections to read 

(^w pq (k)(iV<t uvi)p iVf UV i )q + iV^^iV^J (C/Vt)(t) ® . . . ® (UV^)_JE) (98) 



where the overall sign and normalization is fixed by real virtual cancellation: only with this choice do they 
add up to an infrared finite Fokker-Planck Hamiltonian, which happens to coincide with the expression 
conjectured in Sec. [21 Indeed, the operator appearing in the sum of l|97|) and 198|) is nothing but —H ng of 
Eq. (|34afl with the /W completely determined. In particular, there is no room for any N c dependence of the 
coefficients from this argument. 

To summarize, it has been shown that 

Ed E (e at ^ uv ' +at < uv ^^) uv ^E) = -( H ag ^ uv^ +i2tr(uv^ j } (E) . (99) 
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The only argument still missing is about how to reverse the step leading from l|18|) to 12111 . 

This follows from the fact that the identity holds for arbitrary j and j\ i.e. the fact that the above is 
completely independent of the type of correlator considered. The result is, as advertised in Sec. El the 
evolution equation 

Ed E Z E [UV^} = -H ng Z E [UV^] (100) 

for the weight Z E [UV>] used to define the averages (. . .) u y,{E). This completes the derivation of Eq. (|36|1 . 

I will close this section with the explicit expressions for the non-factorized version of the BMS equation as 
it emerges from this discussion, to point out how the N c limit comes about and to illustrate once more the 
real-virtual cancellations with a concrete example. As this is concerned with emissions from a bare qq jet, 
one needs the operator and probabilities 



G pq :=tv((UV%(UV%)/N c (101a) 
1 S_ __6_ 

N c iS(j p )ij iS(j p )ij 



G pq (E) :=(G pq } UV i(E) = - - \N E \j\j] . (101b) 



The resulting equation then is 

Ed E {tx{{UV%{UV\)/N c ) uv ,{E) 

— -w pq (k) {u(k)(UV% 2 2C t - )uvt ■ (1 ° 2) 

The first term on the r.h.s., carrying the factor u(fc), originates from the real contribution, the second from 
the virtual one. By l|29(l this is equivalent to 

Ed E G pq (E) = J ^^. Wpq (k)(u(k)(G pk G kq N c - ^p) - 2C i G pq ){E) (103a) 

^a s w pq {k)(u{k)(G pk G kq - G pq ) - ^(1 - u(k))G pq )(E) , (103b) 

which in the large N c limit factorizes into the analogues of the BK that triggered these explorations, the BMS 
equation Eqns. © and respectively. To see this explicitly, use 1)84(1 to represent the averaging involved in 

G pq {E) = {G pq ) uv t (E) = (A$ q [U, mj q [U, S]t) H (104) 

and then repeat the argument of Eq. (50). It is worth emphasizing that the real- virtual cancellation 
argument of BMS restated in Sec. applies already to Eq. (|103H itself, without any reference to the fac- 
torization argument. Infrared finitcness in the Sudakov term follows directly from the argument of BMS 
while in the other term it is the definition of G that saves the day. This provides an illustration of how the 
more general reasoning above will lead to infrared finite results via real-virtual cancellation in any evolution 
equation for color singlet objects contained in the new evolution equation. 



8 Conclusions 

The main result of this paper is the evolution equation l(36|l that generalizes the BMS equation to finite N c . 
The result is functional and represents an infinite hierarchy of coupled equations just as its small x analogue, 
the JIMWLK equation. As in this case, a Langcvin formulation has been derived that hopefully can be 
implemented efficiently as a numerical simulation. 

More subtle benefits lie in the method used to prove the result. This was based on a formulation of strongly 
ordered soft emission amplitudes using shower operators that can be used to generate soft gluon clouds 
around any hard seed. The ordered nature of these is the main calculational benefit that eventually allowed 
the derivation of the new evolution equation. This concept should prove useful in most of the new applications 
briefly touched upon below. 
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The amplitudes were then translated into probabilities and several forms for taking the expectation values 
were given. Each of these has its own advantage as shown by their use in the derivation of the evolution 
equation. This in the end has allowed to derive this equation by simply differentiating the expressions found 
with surprisingly simple manipulations. The simplicity of these steps allows the modification of the equation 
to applications with other global and non-global jet observables. 

By the close analogy with JIMWLK and BK, some of the results obtained in the asymptotic regime, such 
as the presence of the saturation scale clearly seem to map on universal features of the jet evolution, such as 
the existence of # cr it discussed by BMS. Further explorations along these lines appear to be promising. 

Conversely, it should be possible to use the strategy employed here to rederive the JIMWLK equation. In 
this case, quantities that are less inclusive than the dipole cross section in DIS should become accessible. 
Diffraction or the formulation of observables that focus in on small impact parameters would emerge by 
formulating appropriate phase space constraints. The definition of observables that are infrared safe regardless 
of impact parameter and gluon density encountered is another item on this list. Such a step will certainly 
come at the price of an additional term in the evolution equation just like the Sudakov term in BMS, but 
would put applications of that equation on solid ground even for small targets like the proton. This should 
allow a new level of insight for the application of saturation ideas to HERA data and ultimately also help 
with applications to RHIC and LHC physics. 

Acknowledgements: I am grateful to Andreas Schafer and Vladimir Braun for a thorough discussion 
of content and presentation and want to thank Einan Gardi, Andreas Freund, Mrinal Dasgupta and Kari 
Rummukainen for a string of communications and discussions on various subtopics. I also want to acknowledge 
partial support by BMBF. 

A Color algebra of the induction step 

The color algebra pattern can be best understood by looking at the small values for n. This uncovers the 
pattern for the general case. Notation below is simplified in order to expose the structure of terms: 

• n=0: starting from U^U q one arrives at 

JpaUpt Uq + JaqUpt Uq = Upt UqApqiJt) 

• n=l: starting from (U^t ai U q )U aibl one finds 

~J pa Uj,t a t ai UqU aibl + J aia {Ul[t a ,t^]U q )U a ^ + J qa Ult a H a U q )U a ^ 
= (J ai a - J pa )Ult a t^UqU a ^ + {J qa ~ J aia )Ult a H a UqU a ^ 

• n=2: starting from (U^t a H a2 U q )U aibl U a2b2 one obtains 

- J pa {Ult a t a H a2 UqU aibl 

+ J ai a(Ul,[t a ,t ai }t a2 Uq)U aibl + Ja 2 a{U^[t a ,t a2 }U q )U aibl 

+ J qa (u^ a H a n a u q ))u aibi 

The general step is now obvious: to collect terms with the same color structure, the first term (from the 
U') pairs with the first term of the first commutator, then iteratively the second term of one commutator 
with the first term of the next. This continues all the way through until the last remaining term of the last 
commutator pairs up with the term of the U factor. In all the pairings, the signs will alternate, yielding the 
Aw(k). Setting the labels accordingly then leads to the induction step given in Eq. (|7T)l . 
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